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The  aim  of  the  research  on  the  Delft  jet-in-hot  coflow  (DJHC)  burner  is  to  gain  understanding  in  the  inter¬ 
play  of  turbulence  and  chemistry  in  conditions  as  encountered  in  devices  operating  in  flameless  combus¬ 
tion  mode,  and  to  test  the  validity  of  numerical  models  when  applied  to  these  flameless  combustion 
conditions.  Datasets  on  velocities,  temperatures  and  qualitative  OH  data  of  several  Dutch  natural  gas 
flames  in  the  DJHC  burner  have  been  obtained  and  are  discussed  in  this  paper.  It  was  found  that  the  mean 
velocity  and  turbulent  stresses  are  not  significantly  affected  by  the  chemical  reactions,  which  is  in  line 
with  the  very  moderate  increase  of  mean  temperatures  in  the  flames.  Even  at  heights  where  flame  struc¬ 
tures  are  present,  peak  temperatures  do  not  always  approach  the  adiabatic  flame  temperature.  With  both 
flame  luminescence  and  OH-PLIF  measurements,  it  is  seen  that  chemical  reactions  begin  to  occur  at  a 
lower  location  when  the  jet  velocity  (and  thereby  the  jet  Reynolds  number)  is  increased.  By  analysing 
the  velocity  and  temperature  data  in  the  near-nozzle  region,  the  entrainment  of  coflow  fluid  into  the  tur¬ 
bulent  jet  has  been  quantified.  The  increased  entrainment  of  the  higher  Reynolds  number  jet,  in  combi¬ 
nation  with  the  positive  temperature  gradient  in  radial  direction  in  the  near  field  of  the  jet,  is  shown  to  be 
responsible  for  the  decrease  of  the  height  where  reactions  start  to  occur. 
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1.  Introduction 

By  operating  furnaces  and  boilers  [1  j  and  also  gas  turbines  [2]  in 
flameless  (or  MILD)  combustion  mode,  high  efficiencies  can  be 
combined  with  very  low  NOx  emissions.  Flameless  combustion 
therefore  offers  attractive  characteristics,  from  both  an  environ¬ 
mental  and  an  economic  point  of  view,  even  more  so  in  the  light 
of  the  increasingly  strict  regulations  on  C02  and  NOx  [3]  emissions. 
However,  large  scale  implementation  of  this  technique  is  ham¬ 
pered  by  the  risks  associated  with  adopting  a  relatively  unknown 
new  technique. 

The  ability  to  predict  (turbulent)  flames  by  computer  models 
can  speed  up  the  development  and  implementation  of  more  ad¬ 
vanced  combustion  techniques  [4].  A  general  insight  in  the  specific 
demands  of  numerical  models  to  accurately  predict  the  character¬ 
istics  (for  instance,  heat  fluxes  and  emissions)  of  flameless  com¬ 
bustion  would  therefore  be  beneficial  to  its  implementation. 
Several  computational  studies  have  been  carried  out  on  (semi)- 
industrial-scale  setups  operating  in  flameless  mode  [5-8].  How¬ 
ever,  dealing  with  the  complex  turbulence-chemistry  interaction 
in  this  combustion  mode  accurately  and  computationally  effi¬ 
ciently  remains  a  challenge.  Moreover,  the  reasons  for  success  or 
failure  of  particular  models  are  still  unclear.  This  is,  at  least 
partially,  caused  by  the  limited  (optical)  access  provided  by  these 
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setups,  hindering  detailed  studies  on  for  instance  the  statistics  of 
the  turbulent  flow  and  temperature  field  throughout  the  domain. 
As  a  result,  there  are  not  many  quantities  to  be  used  for  validation. 
Given  the  large  number  of  physical  processes  occuring  in  a  furnace 
environment,  this  makes  the  isolation  and  investigation  of  a  single 
aspect  troublesome. 

To  study  the  interplay  of  turbulence  and  chemistry  in  more  de¬ 
tail,  it  is  therefore  advantageous  to  simplify  the  geometry,  while 
retaining  the  physical  processes  of  interest,  and  to  make  the  setup 
optically  accessible.  Research  on  such  prototypical  laboratory-scale 
flames  with  the  level  of  detail  offered  by  modern  laser-diagnostics 
has  resulted  in  progress  in  the  area  of  turbulent  combustion,  by 
providing  data  for  model-validation  and/or  by  the  information 
contained  in  the  experimental  data  itself  [4,9]. 

In  this  light,  studying  flameless  combustion  in  an  open,  uncon¬ 
fined  setup  might  give  valuable  insights,  provided  that  the  circum¬ 
stances  resemble  those  found  locally  in  a  furnace.  Several  laboratory 
scale  jet  flames  with  a  vitiated  coflow  have  been  developed,  most 
notably  the  Cabra  et  al.  burner  [10],  and  the  Adelaide  burner  [11], 

The  Cabra  burner  was  designed  to  capture  the  features  of  recir¬ 
culation  burners.  Its  most  striking  feature  is  the  large  diameter  of 
its  coflow,  providing  homogeneous  boundary  conditions  over  a 
large  range,  such  that  it  can  be  considered  a  two-stream  problem. 
This  burner  was  not  aimed  specifically  at  investigating  flameless 
combustion,  and  it  might  be  argued  that  the  oxygen  mole  fraction 
in  the  coflow  (X02  of  0.15  or  0.12,  for  the  H2-N2  and  CH4/air  jet 
flames  respectively  [12])  is  rather  high  for  that  purpose.  Multiple 
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detailed  measurement  campaigns  have  been  carried  out,  providing 
a  comprehensive  dataset  of  single-point  scalar  statistics  [13,14] 
and  information  on  velocities  [15],  Furthermore,  joint  temperature, 
OH,  and  H2CO  measurements  have  been  performed  [16], 

The  Adelaide  burner  was  designed  to  study  flameless  combus¬ 
tion,  with  low  oxygen  mass  fractions  in  the  coflow  (typically 
3-9%)  and  a  small  temperature  increase  in  the  reaction  zone  as  a 
result.  Detailed  species  measurements  [11]  and  also  joint  temper¬ 
ature,  OH,  and  H2CO  measurements  [17]  on  flames  with  CH4-H2  in 
the  fuel  jet  have  been  carried  out. 

in  this  paper  the  detailed  results  of  temperature  and  velocity 
measurements  in  the  Delft  jet-in-hot  coflow  burner  are  reported, 
for  three  different  coflow  cases  (different  temperatures  and  oxygen 
mass  fractions)  and  jet  Reynolds  numbers  ranging  from  2500  to 
8800.  The  results  are  analysed  with  the  previous  findings  on  the 
flame  stabilisation  in  mind,  namely,  the  flame  is  stabilised  by  ran¬ 
domly  formed  autoignition  spots  that  grow  while  being  convected 
downstream  by  the  mean  flow  [18].  Datasets  of  the  various  flames, 
containing  information  on  velocities,  temperatures,  OH-fields  and 
the  coflow-composition  are  available  upon  request. 

2.  Experimental  setup 

2.1.  DJHC-burner 

The  design  of  the  Delft  jet-in-hot  coflow  (DJHC)  burner  is  based 
on  that  of  the  Adelaide  jet-in-hot  coflow  burner.  The  main  differ¬ 
ence  is  that  the  latter  uses  addition  of  N2  to  cool  down  the  coflow 
whereas  the  DJHC  burner  uses  cooling  of  the  coflow  through  radi¬ 
ative  and  convective  heat  losses  along  the  burner  pipe.  Another  dif¬ 
ference  is  found  in  the  design  of  the  secondary  burner,  which 
allows  for  the  addition  of  seeding  particles  in  the  DJHC-burner, 
to  act  as  tracers  in  for  instance  LDA  or  P1V  measurements.  The 
schematic  of  the  burner  is  shown  in  Fig.  1 .  The  secondary  burner, 
which  generates  the  hot  coflow,  is  a  ring-burner.  Its  design  was 
found  to  yield  an  axi-symmetric  flame,  and  it  is  operated  in  par¬ 
tially  premixed  mode  (by  adding  24  nl/min  of  air  to  the  fuel 


Fig.  1.  Schematic  design  of  the  Delft  JHC  burner,  with  relevant  dimensions  in  the 
inset.  Two  dimensions  are  given  for  the  vertical  extent  of  the  fuel  tube  above 
the  outer  tube,  for  the  cold  situation  and  the  hot  situation.  The  difference  is 
due  to  the  uneven  thermal  expansion  of  the  fuel-  and  outer  tube. 


stream)  to  help  it  stabilise.  The  secondary  burner  creates  a  flame 
with  a  length  of  about  0.4  m.  A  grid  is  located  0.11  m  upstream 
of  the  burner  exit  plane,  to  keep  the  fuel  pipe  centered  and  to  help 
the  coflow  lose  heat.  The  outer  burner  tube  radiates  most  strongly 
at  the  height  of  the  distribution  grid,  indicating  that  this  cooling 
mechanism  is  effective.  The  central  fuel  pipe  (with  internal  diam¬ 
eter  4.5  mm)  is  cooled  by  constantly  flushing  air  through  the  con¬ 
centric  cooling  air  ducts,  thus  preventing  excessive  heating  of  the 
main  fuel  jet.  The  main  flow  of  air  passes  through  the  air  inlets 
at  the  bottom  and  it  is  this  flow  that  carries  the  seeding  particles. 
A  z,  r-coordinate  system  is  used  centered  on  the  exit  of  the  fuel 
pipe.  The  velocities  corresponding  to  these  directions  are  u  and  v, 
respectively. 

2.2.  LDA  system 

LDA  measurements  were  performed  with  a  two-component, 
dual  beam  TSI-system.  The  green  line  (514.5  nm)  and  blue  line 
(488  nm)  of  a  10  W  Continuum  Argon-ion  laser  were  used  to  mea¬ 
sure  the  axial  and  radial  velocity  components  directly.  Two  of  the 
incident  beams  (one  of  each  color)  were  frequency  pre-shifted  over 
40  MHz  by  a  Bragg  cell  to  enable  the  detection  of  instantaneous 
flow  reversals  and  stagnant  flow.  The  focusing  lens  had  a  82  mm 
aperture  and  a  focal  length  of  250  mm.  The  length  and  diameter 
of  the  measurement  volume  were  1.7  mm  and  0.12  mm,  respec¬ 
tively.  The  fringe  distances  for  the  green  and  blue  channel  were 
2.6  pm  and  2.5  pm.  Alumina  (A1203)  particles  with  an  average  size 
of  about  1  pm  were  used  as  seeding  particles.  Two  cyclone-type 
particle  generators  were  used  to  seed  the  air  and  fuel  separately. 
The  generators  have  a  provision  that  enables  the  control  of  the 
seed  density  in  both  the  coflow  and  in  the  fuel  jet.  This  provision 
was  used  to  equalise  the  seeding  density  in  the  both  flows,  thereby 
minimizing  errors  related  to  so-called  conditional  seeding.  The 
data  rate  (number  of  bursts  per  second)  was  used  as  an  indicator, 
as  it  is  proportional  to  the  fluid  density,  the  velocity  magnitude  and 
the  seeding  density.  This  proportionality  was  therefore  used  to 
equalize  the  seeding  densities:  the  ratio  of  the  data  rates  of  the  un¬ 
mixed  fuel-  and  coflow  streams  was  made  approximately  equal  to 
the  ratio  of  their  products  of  axial  velocity  and  fluid  density.  The 
light  scattered  by  the  seeding  particles  was  collected  in  back-scat- 
ter  mode.  The  photomultiplier  output  signals  were  electronically 
down-mixed,  and  subsequently  fed  to  a  FSA-3000  signal  processor 
to  determine  the  instantaneous  velocity  of  light-scattering  parti¬ 
cles.  All  statistics  were  computed  as  transit-time  weighted  results 
to  eliminate  the  effects  of  the  velocity  bias.  Autocorrelation  func¬ 
tions  of  the  axial  velocity  component  were  constructed  from  time 
series  with  4  x  105  velocity  samples  that  were  acquired  at  a  mean 
data  rate  of  approximately  500  Hz  by  using  the  slotting  method 
with  local  normalisation  [19]. 

2.3.  CARS  system 

Temperatures  were  determined  with  a  CARS  system  that  has 
been  described  in  detail  elsewhere  [20],  It  is  based  on  an  injec¬ 
tion-seeded,  frequency-doubled  Nd:YAG  laser  (Spectron  SL805 
SLM),  which  yields  500  mj  per  pulse  at  532  nm  with  a  pulse  dura¬ 
tion  of  12  ns  at  10  Hz  repetition  rate.  About  80%  of  the  radiation  is 
used  to  pump  a  modeless  Stokes  dye  laser  (Mode-X  ML-3),  emit¬ 
ting  a  broadband  profile  around  607  nm  for  Rhodamine  640  in 
methanol.  The  remaining  20%  of  the  pump  laser  travels  along  a  de¬ 
lay  line  and  is  split  into  two  beams  with  equal  intensity.  The  Stokes 
beam  and  the  two  beams  at  532  nm  are  focused  by  an  aplanat  lens 
with  a  focal  length  of  300  mm  in  a  planar-boxcars  phase-matching 
configuration.  With  this  configuration,  a  CARS  probe  volume  of 
700  pm  length  and  35  pm  diameter  is  obtained.  The  generated 
CARS  radiation  is  recollimated  and  combined  with  an  attenuated 
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sample  of  the  Stokes  beam  on  a  dichroic  beam  splitter.  The  beams 
are  focused  onto  the  entrance  of  an  echelle  spectrometer,  and  their 
spectra  are  dispersed  on  a  CCD  detector  with  1100  x  330  pixels. 
The  spectra  are  contained  in  two  strips  of  1100  intensity  values, 
which  are  digitized  by  an  18-bit  AD  converter  and  stored.  With 
Dacapo  software  [20]  the  CARS  spectrum  is  referenced  to  the 
simultaneously  measured  Stokes  excitation  profile  and  fitted  to  a 
library  of  theoretical,  temperature-dependent  spectra.  The  single¬ 
shot  imprecision  of  the  system  is  1  -4%  over  a  range  from  2000  K 
to  300  K.  The  inaccuracy  is  estimated  to  be  20  K.  For  each  point 
in  space,  mean  temperatures  were  determined  from  the  results 
of  1000  single-shot  CARS  spectra. 

2.4.  OH-PLIF  imaging  system 

A  frequency-doubled,  pulsed  Nd:YAG  laser  (Spectra  Physics 
PR0250-10)  pumped  a  Syrah  PrecisionScan  dye  laser  using  Rhoda- 
mine  590  in  methanol.  The  second  harmonic  radiation  from  the 
dye  laser  was  used  to  excite  the  Qi(6)  line  of  the  A2X+-X2n 
(1,0)  transition  at  282.927  nm.  The  UV  energy  was  measured  to 
be  about  10  mj/pulse.  Three  cylindrical  lenses  transformed  the  la¬ 
ser  beam  into  a  sheet  with  a  height  of  approximately  80  mm  and 
an  approximate  thickness  of  200  pm.  Fluorescence  from  the 
A2£+  -  X2n  (0, 0)  and  1,1)  transitions  at  305-3 1 5  nm  was  collected 
through  a  Semrock  narrowband  filter  and  a  UV-Nikkor  105  mm 
// 4.5  lens.  The  images  were  recorded  by  an  intensified  high-speed 
camera  (Lambert  Instruments  HI-CAM  CR)  operating  at  the  full  res¬ 
olution  of  1280  x  1024  pixels.  A  mean  off-resonance  image  was 
subtracted  from  the  images  to  determine  the  mean  fluorescence 
signal.  Given  a  certain  point  in  the  flame,  the  mean  OH  fluores¬ 
cence  signal  depends  on  the  laser  sheet  intensity  and  the  sensitiv¬ 
ity  of  the  camera.  This  aspect  was  investigated  by  traversing  the 
flame  vertically  through  the  light-sheet  keeping  the  camera  posi¬ 
tion  constant  and  determining  the  mean  signal  strength.  It  was 
found  that  the  resulting  mean  signal  strength  for  a  given  point  in 
the  flame  (flame  DJHC-I,  at  z  rj  90  mm)  is  within  30%  of  the  value 
in  the  center  of  the  light  sheet  over  a  height  of  55  mm,  which  is 
considered  acceptable  for  the  purpose  of  flame  visualisation.  To  ob¬ 
tain  statistics  of  sufficient  quality,  1000  images  were  recorded  for 
each  flame. 

2.5.  Flue  gas  measurements 

Oxygen  measurements  were  performed  with  a  Testo  335  flue¬ 
gas  analyser,  with  a  specified  inaccuracy  of  ±0.20%.  The  measured 
oxygen  volume  fractions  in  the  coflow  were  converted  to  mass 
fractions  using  the  results  from  equilibrium  chemistry  calculations 
with  the  species  of  the  detailed  Warnatz-mechanism  that  includes 
C-2  chemistry  (see  [21],  pp.  333-339). 

3.  Case  description 

Three  different  settings  of  the  secondary  burner,  affecting 
mainly  the  coflow  temperature,  were  used.  The  resulting  coflow 


properties  are  summarised  in  Table  1.  The  coflow  properties 
(velocity,  temperature  and  oxygen  mass  fraction)  vary  over  the 
cross-section.  The  table  gives  the  maximum  of  the  mean  tempera¬ 
ture  profile  measured  at  z  =  3  mm.  The  reported  02  mass  fraction  is 
mass  flux-averaged. 

Under  adiabatic  and  isobaric  circumstances  the  total  (sensible 
plus  chemical)  enthalpy  is  unaffected  by  chemical  reaction,  and 
as  such  it  is  a  conserved  quantity.  The  heat  losses  of  the  coflow 
bring  about  an  enthalpy  deficit:  a  difference  between  the  adiabatic 
total  enthalpy  and  the  actual  total  enthalpy.  As  the  enthalpy  deficit 
of  the  coflow  is  carried  to  the  flame  zone,  lower-than-adiabatic 
flame  temperatures  (of  cold  air  with  cold  fuel)  are  expected.  The 
enthalpy  deficit  of  a  stoichiometric  mixture  of  fuel  and  an  oxidiser 
that  consists  of  flue  gases  can  be  expressed  as  a  function  of  the 
mixture  fraction  of  the  flue  gases  %°,  the  stoichiometric  mixture 
fraction  £st  (both  mixture  fractions  based  on  the  fresh  fuel  and 
air  streams),  and  the  enthalpy  deficit  of  the  oxidiser  Afi°: 

mst)  =  \^Ah°.  (i) 

i  -  r 

Because  more  heat  is  lost  by  radiation  of  the  outer  tube  when  the 
coflow  is  hotter,  the  adiabatic  flame  temperature  is  lowest  in  the 
DJHC-I  case.  The  adiabatic  flame  temperatures  reported  in  Table  1 
are  based  on  equilibrium  chemistry  calculations,  using  the  mea¬ 
sured  coflow  properties  at  r  =  20  mm. 

The  fuel  used  in  the  central  jet  was  Dutch  natural  gas  which  has 
an  approximate  composition  of  81%  methane,  4%  ethane,  14% 
nitrogen  and  1%  higher  alkanes,  by  volume.  Its  adiabatic  flame 
temperature  in  ambient  air  is  2210  K.  Different  mass  flows  were 
used  in  the  fuel  jet,  namely  10.7,  16.1  and  30.0nl/min  resulting 
in  jet  Reynolds  numbers  between  2500  and  8800,  based  on  the 
bulk  velocity  of  the  fuel  and  the  density  and  dynamic  viscosity  of 
Dutch  natural  gas  at  the  estimated  fuel  temperature  at  the  nozzle 
exit.  These  temperatures  have  been  determined  using  mass  conser¬ 
vation,  by  comparing  the  measured  mass  flow  and  the  integrated 
volume  flow  at  the  nozzle  exit,  from  which  the  density  and  thus 
the  temperature  follows.  The  reason  for  taking  this  approach  rather 
than  using  the  temperature  data  directly  is  the  quality  of  the  CARS 
signal  obtained  close  to  the  fuel  pipe  exit.  Here,  the  N2-CARS  spec¬ 
trum  is  obscured  by  the  strong  non-resonant  methane  signal.  In 
addition,  vibrational  CARS  suffers  from  increasing  imprecision  at 
lower  temperatures. 

The  jet  Reynolds  numbers  are  somewhat  limited  in  magnitude. 
Though  this  is  generally  not  advantageous  for  RANS  calculations,  it 
might  on  the  other  hand  be  attractive  for  LES  studies.  Measure¬ 
ments  were  done  on  a  selection  of  the  possible  combinations  of 
coflow  cases  with  jet  mass  flows.  The  studied  combinations  are 
indicated  in  Table  2,  along  with  their  jet  Reynolds  numbers  and 
fuel  temperatures.  The  higher  and  lower  values  of  the  jet  Reynolds 
number  were  only  fully  studied  with  LDA  and  CARS  in  the  DJHC-I 
flames. 

The  profiles  of  mean  axial  velocity  of  the  hot  coflow  and  the  fuel 
jet  are  shown  in  Fig.  2.  The  profiles  of  temperature  and  oxygen  con¬ 
tent  are  shown  in  Figs.  3  and  4.  The  jet  Reynolds  number  does  not 


Table  1 

Characteristics  of  the  coflow  stream  of  the  studied  DJHC  flames.  The  amount  of  natural  gas  to  generate  the  coflow  and  the  total  amount  of  air  (premixed  and  non-premixed)  are 
given  in  columns  two  and  three.  rmax;co  is  the  maximum  coflow  temperature,  whereas  Vo2:av  is  a  mass  flux  weighted  average  between  r  =  2.5  mm  and  r  =  35  mm,  at  z  =  3  mm.  The 
subscript  “20”  denotes  properties  evaluated  at  r  =  20  mm.  The  mixture  fraction  is  based  on  the  inlet  streams  of  the  burner  (normal  air  and  natural  gas),  whereas  f[t;20  is  the 
rescaled  stoichiometric  mixture  fraction,  of  a  mixture  with  the  coflow  at  r  =  20  mm  as  oxidiser  and  natural  gas  as  fuel.  The  uncertainty  of  the  last  digit  in  the  values  for  enthalpy 
loss  and  mixture  fraction  are  given  in  parentheses. 

Sec.  fuel  Tot.  air  Tmax  y02;av  A h°0  Gd;2o 

_ (nl/min) _ (nl/min) _ (K) _ (-)(%) _ (J/kg) _ (-) _ H _ (K) _ 

DJHC-I  16.1  224  1  540  7.6  -4.4(5)  x  105  0.052(1)  0.019  1  950 

DJHC-V  15.3  231  1460  8.8  -4.0(5)  x  105  0.048(1)  0.024  1  980 

DJHC-X  14.2  239  1395  10.9  -2.8(5)  x  105  0.042(1)  0.029  2050 
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Table  2 

Jet  Reynolds  numbers  with  fuel  temperatures  between  parentheses,  of  the  combi¬ 
nations  of  different  coflows  and  jet  mass  flows  that  have  been  fully  studied,  with  both 
LDA  and  CARS.  For  Dutch  natural  gas,  one  normal  liter  corresponds  to 
0.833  x  10_3kg.  The  relative  uncertainty  in  the  fuel  temperature  is  estimated  at 
±5%,  the  uncertainty  in  the  Reynolds  numbers  is  estimated  at  ±200  for  the  Re  =  2500 
and  4100  cases,  and  ±300 1<  for  the  other  cases.  The  default  traverses  are  at  z  =  3, 15, 
30,  60,  90,  120  and  150  mm  (radial  traverses)  and  a  centerline  traverse.  The  cases 
superscripted  with  1  have  a  reduced  number  of  measurement  locations  with  the 
traverses  at  z=  15,  30  and  150  mm  and  the  centerline  traverse  omitted. 


Mjet  (nl/min) 

10.7 

16.1 

30.0 

DJHC-I 

2500  (470  K)1 

4100 (430  K) 

8800  (360  K) 

DJHC-V 

4600  (380  K)1 

DJHC-X 

4600 (380  K) 

influence  the  oxygen  mass  fraction  nor  the  temperature  signifi¬ 
cantly  at  z  =  3  mm.  In  both  the  temperature  and  oxygen  mass  frac¬ 
tion  profiles,  asymmetries  are  seen  that  exceed  the  inaccuracy  of 
the  measurements.  The  maximum  difference  between  the  oxygen 
mass  fraction  at  the  left  and  right  side  is  0.8%,  whereas  the  left- 
right  deviation  of  the  peak  temperature  is  between  20  and  40  K. 
With  the  current  design,  asymmetries  of  this  magnitude  cannot 
be  avoided. 


Fig.  2.  Axial  velocities  at  z  =  3  mm,  as  measured  with  LDA.  The  inset  shows  the 
turbulence  intensity,  calculated  as  /  =  (u'u1  +  2z/z/)1/2/U0,  with  U0  the  centerline 
velocity  at  the  nozzle  exit. 


Fig.  3.  Mean  and  RMS  of  temperature  at  z  =  3  mm,  measured  with  CARS. 


4.  Visual  and  OH-LIF  observations 

The  visual  appearance  of  several  DJHC  flames  is  shown  in  Fig.  5. 
As  a  reference,  the  bottom  right  image  shows  a  conventional  lifted 
flame  with  a  jet  Reynolds  number  of  5000  and  a  coflow  velocity  of 
0.5  m/s.  All  pictures  were  taken  with  identical  exposure  time,  CCD 
sensitivity  and  aperture.  The  jet-in-hot-coflow  flames  have  consid¬ 
erably  weaker  flame  luminescence.  In  the  DJHC-X  flame,  hardly 


Fig.  4.  Oxygen  mass  fractions  at  z  =  3  mm,  measured  with  a  flue  gas  analyser. 


(a)  DJHC-I  Re  2.5K  (b)  DJHC-I  Re  4.  IK  (c)  DJHC-I  Re  8.8K 


(d)  DJHC-V  Re  4.6K  (e)  DJHC-X  Re  4.6K  (f)  lifted  Re  5K 


Fig.  5.  Images  of  the  DJHC-flames  and  a  conventional  lifted  flame,  all  with  0.5  s 
exposure  time  and  identical  aperture.  The  top  of  the  photograph  is  in  each  case  at 
around  z=  180  mm.  The  coordinate  system  is  indicated  in  Fig.  5a. 


E.  Oldenhof  et  al. / Combustion  and  Flame  158  (201 1 )  1553-1563 


1557 


Fig.  6.  Lift-off  heights  (continuous  lines)  and  location  of  first  occurrence  of  ignition 
kernels  (dotted  lines)  as  a  function  of  the  jet  mass  flow  (in  nl/min).  Note  that  the 
flame  DJHC-X  has  a  lift-off  height  exceeding  150  mm  in  all  cases. 

any  flame  luminescence  is  seen  below  z  =  1 50  mm  suggesting  that 
little  chemical  reactions  occur  in  this  region.  This  is  confirmed  by 
OH-PLIF  observations,  and  it  makes  this  flame  suitable  as  a  com¬ 
parison  case  to  the  reacting  flames. 

Figure  6  shows  the  lift-off  heights  based  on  chemiluminescence 
images  obtained  with  the  intensified  high-speed  camera,  with  the 
procedure  outlined  in  [18].  The  lift-off  heights  initially  show  a 
decreasing  trend  as  function  of  the  jet  Reynolds  number.  A  similar 
trend  can  be  seen  in  the  Adelaide  JHC  burner  in  the  cases  with  3% 
oxygen  in  the  coflow  [17].  It  has  been  demonstrated  that  the  de¬ 
crease  in  lift-off  height  is  related  to  the  lower  location  of  first  occu¬ 
rence  of  ignition  kernels,  denoted  with  zb;min  [18].  Both  the 
variation  in  lift-off  height  hi  and  zb:min  with  the  jet  Reynolds  num¬ 
ber  are  shown  in  Fig.  6.  A  striking  fact  is  that  the  initial  decrease  for 
all  cases  is  identical.  Apparently,  the  mechanism  responsible  for 
lowering  the  lift-off  heights  has  equal  strength  in  all  three  coflow 
cases.  It  is  therefore  expected  that  changes  in  the  flow  field  are 
responsible.  This  aspect  will  be  investigated  further  in  Section  7. 

The  OH-PLIF  images  in  Fig.  7  were  constructed  by  taking  the 
root  mean  square  (RMS)  of  the  recorded  intensities.  By  taking  the 
RMS  value  of  the  OH-PLIF  signal,  rather  than  the  mean,  the  occa¬ 
sional  appearance  of  flame  pockets  at  the  base  of  the  flame  is 
clearly  highlighted.  The  axial  location  where  flame  structures  begin 
to  affect  significantly  the  RMS  value  of  the  OH  signal,  decreases  by 
approximately  20  mm  from  Rejet  =  2500  to  Rejet  =  8800  in  flame 
DJHC-I.  This  is  consistent  with  the  earlier  observations  based  on 


Fig.  8.  Mean  axial  velocities  at  z  =  3,  60  and  120  mm  for  the  three  different  coflow 
cases,  at  a  jet  Reynolds  number  of  4100  (DJHC-I)  and  4600  (DJHC-V  and  DJHC-X). 

flame  luminescence  regarding  the  influence  of  the  jet  Reynolds 
number  on  zb  min.  The  mean  OH  signal  at  z  =  90  mm  peaks  at  a  ra¬ 
dius  of  approximately  13  mm  for  the  lowest  jet  Reynolds  number 
and  at  approximately  12  mm  for  the  highest  jet  Reynolds  number. 

5.  Flow  field 

5.1.  Mean  velocities  and  Reynolds  stresses 

Figure  8  shows  the  mean  axial  velocity  profiles  for  the  three  co¬ 
flow  cases  with  a  jet  mass  flow  of  16.1  nl/min  at  z  =  3  mm,  60  mm 
and  120  mm  (z/d  =  2/3,  13  1/3  and  26  2/3,  respectively).  Although 
flame  DJHC-1  has  a  lower  lift-off  height  than  flame  DJHC-X  (ob¬ 
served  in  both  CARS  measurements  and  visually),  this  is  not  re¬ 
flected  in  a  difference  in  the  mean  velocity  profiles.  From  this  it 
can  be  concluded  that  the  effect  of  combustion  on  the  mean  flow- 
field  in  the  lower  portion  of  the  flame  is  minimal.  In  Fig.  9,  the  cen¬ 
terline  decay  of  the  mean  axial  velocity  profile  (normalised  by  the 
jet  exit  velocity  U0  for  two  coflow  cases  and  different  jet  Reynolds 
numbers)  are  presented.  Whereas  the  difference  in  velocity  field 
between  the  coflow  cases  is  small,  there  is  a  clear  distinction  be¬ 
tween  the  low-  and  high  jet  Reynolds  number  case.  The  low 
Reynolds  number  case  shows  a  weaker  initial  decay  but  catches 
up  later,  which  is  also  seen  in  the  normalised  value  of  the  Reynolds 
normal  stress  u'u'. 

The  location  of  maximum  shear  stress  u'v'  is  a  useful  indicator 
of  the  jet  width.  In  the  self-similar  region  of  a  jet,  it  peaks  at  a 
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Fig.  7.  RMS  values  of  the  OH-fluorescence  signal  in  the  flame  stabilisation  region,  DJHC-1,  ReJet  =  2500  (a),  ReJet  =  4100  (b)  and  ReJet  =  8800  (c). 
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Fig.  9.  Centerline  decay  of  the  mean  axial  velocity,  all  cases. 


x  10  3 


Fig.  10.  The  shear  stress  u'v'  at  radial  traverses  at  z  =  60  and  120  mm,  all  cases. 


radius  of  roughly  0.7  times  the  jet  half-width  (the  radius  where  the 
mean  axial  velocity  is  half  that  at  the  centerline)  [22].  Atz  =  60  mm 
and  120  mm,  the  width  of  the  jet  of  different  cases  is  very  compa¬ 
rable,  and  the  shear  stress  u'v'  is  of  similar  magnitude  when  nor¬ 
malised  by  the  jet  exit  velocity  squared  (Fig.  10). 

5.2.  Turbulence  length  scales 

As  turbulence  dominates  the  transport  of  scalars,  an  adequate 
description  of  the  turbulent  flow  field  is  important.  The  Eulerian 
velocity  time  scales  can  be  determined  from  a  stationary  LDA 
probe  by  constructing  autocorrelation  functions  of  the  velocity 
components.  The  turbulence  intensity  on  the  jet  centerline  is 
sufficiently  low  (Vu^/U  «  15%)  to  convert  these  time  scales  to 
length  scales  using  Taylor’s  approximation,  as  the  correction 


proposed  by  Lumley  [23]  on  the  Taylor  length  scale  is  smaller 
than  10%.  The  single  point  LDA  measurements  can  be  used  to 
determine,  for  instance,  integral-  and  Taylor  length  scales.  The 
Reynolds  number  based  on  the  longitudinal  Taylor  length  scale 
Xf,  ReT  =  ^u^Xf/v  (with  v  the  kinematic  viscosity)  forms  a  useful 
quantity  to  compare  experimental  and  numerical  data  because  it 
can  be  determined  from  numerical  data  as  well  (LES  or  DNS).  The 
longitudinal  Taylor  length  scale  follows  from  the  following  expan¬ 
sion  of  the  spatial  autocorrelation  function,  as  function  of  the 
axial  separation  Az: 

A/2 

P(Az)  =  1  — +  ©(Az4),  (2) 

Xf 

and  the  transversal  Taylor  length  scale  Xg  is  defined  similarly.  These 
scales  are  related  to  the  mean  square  derivatives 
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To  determine  the  Taylor-scale  Reynolds  number,  the  kinematic 
viscosity  was  determined  based  on  equilibrium  chemistry  calcula¬ 
tions  (on  the  centerline,  the  mixture  is  still  rich  and  the  difference 
between  reacting  and  non-reacting  is  minimal).  The  used 
kinematic  viscosities  at  z=60  and  120  mm  are  8  x  10-5  and 
12  x  1 0-5  m2/s,  respectively.  The  integral  longitudinal  length  scale 
was  determined  by  integration  of  the  longitudinal  autocorrelation 
function. 

The  properties  of  the  flowfield  on  the  centerline  at  z  =  60  and 
z=120mm  are  summarised  in  Table  3.  The  transverse  Taylor 
length  scale  is  smaller  than  the  longitudinal  length  scale.  Note  that 
in  isotropic  turbulence  Xf  =  V2Xg  (see  [22],  p.  199).  The  Taylor 
length  scale  grows  with  increasing  distance  from  the  jet  exit  along 
with  the  integral  length  scale.  The  ratio  between  the  Taylor  and  the 
integral  length  scales  is  not  observed  to  decrease  significantly  with 
the  jet  Reynolds  number.  Scaling  arguments  (see  for  instance  [22], 
p.  200)  and  measurements  (as  reported  for  instance  in  [24])  sug¬ 
gest  an  inverse  square  root  relation,  but  the  difference  in  Reynolds 
numbers  is  too  small  here  to  be  of  significance. 

To  calculate  the  dissipation  of  turbulent  kinetic  energy  in  gen¬ 
eral  non-isotropic  turbulence,  twelve  different  derivative  correla¬ 
tions  need  to  be  measured,  which  is  practically  impossible  with 
single  two-component  LDA.  The  assumption  of  local  axisymmetry 
leads  to  a  relaxation  of  these  demands,  but  derivatives  in  the  radial 
directions  are  still  needed  [25].  As  this  direction  is  not  aligned  with 
the  mean  flow,  these  derivatives  can  still  not  be  accessed  with  a 
stationary  probe.  Therefore,  to  get  an  estimate  of  the  smallest 
scales,  isotropy  of  the  smallest  scales  will  be  assumed.  As  a  mea¬ 
sure  of  the  anisotropy,  Kj  is  also  given  in  Table  3.  This  quantity  is 
defined  as: 


Table  3 

Turbulence  properties  on  the  centerline  at  two  heights  above  the  jet  exit,  determined  with  two-component  LDA  measurements. 
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DJHC-I 

4100 

60 
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5000 
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and  it  is  equal  to  one  in  isotropic  conditions  [25].  As  can  be  seen  in 
the  table,  the  deviations  from  unity  are  not  large  at  z  =  120  mm.  An 
estimate  of  the  dissipation  rate  of  turbulent  kinetic  energy  based  on 
the  expression  for  isotropic  turbulence  (e  =  30 vu'2 /2)  is  included 
in  the  final  column.  The  Kolmogorov  length-scale  XK  =  (v3e~')1/4  is 
in  the  order  of  1  x  10~4  m. 

6.  Temperature  field 

The  high-speed  intensified  imaging  reported  in  [18]  described 
the  probability  of  the  presence  of  flame  structures  as  a  function 
of  axial  height,  for  various  flames.  This  probability  was  expressed 
in  the  parameter  denoted  with  Pbl.  This  parameter  will  be  used 
in  this  section  for  comparison,  a  discussion  on  its  definition  is  pro¬ 
vided  in  the  aforementioned  paper.  The  temperature  measure¬ 
ments  at  z  =  60  mm  (not  shown)  show  no  sign  of  any  reactions. 
On  the  other  hand,  at  z  =  90  mm  a  mild  temperature  increase  can 
be  observed  in  the  DJHC-I  flame,  see  Fig.  11.  This  increase  is  hardly 
noticeable  in  the  mean,  but  is  visible  in  the  temperature  RMS,  and, 
most  notably,  in  the  peak  temperature  Tgg.  This  peak  temperature 
is  defined  as  the  temperature  where  the  cumulative  density  func¬ 
tion  exceeds  99%  (it  is  the  99th  percentile).  Although  flame  struc¬ 
tures  are  present  during  approximately  80%  of  the  time  at 
z  =  90  mm  (Pbl  80%),  the  peak  temperature  Tgg  does  not  exceed 
1770  K,  relatively  far  from  the  estimated  adiabatic  flame  tempera¬ 
ture  of  1950  K.  Comparing  the  locations  of  the  peaks  in  Tgg  in 
Fig.  11  (at  a  radial  location  of  around  12  mm,  which  is  in  agree¬ 
ment  with  the  location  of  the  flame  as  observed  in  the  OH-PLIF 
measurements)  with  the  measured  shear  stress  profile  at 
z  =  90  mm  (which  peaks  at  a  radius  of  5  mm),  the  reaction  zone 
clearly  resides  outside  the  region  of  most  intense  turbulence.  This 
can  be  understood  from  the  value  of  the  stoichiometric  mixture 
fraction,  being  even  lower  than  that  of  a  flame  of  Dutch  natural 
gas  in  air  (around  0.02  instead  of  0.07)  due  to  the  reduced  oxygen 
content  of  the  coflow. 

At  z  =  120  mm,  the  presence  of  the  reaction  zone  of  flame  DJHC- 
1  is  visible  in  the  mean  temperature  profile,  whereas  this  is  less 
clear  in  flame  DJHC-V,  see  Fig.  12.  An  interesting  feature  is  that 
Tgg  of  flame  DJHC-I,  increases  by  some  140  K  from  z  =  90  to 
z  =  120  mm,  to  around  1910  K.  As  the  flame  probability  Pbl  is  al¬ 
ready  close  to  100%  at  z  =  90  mm,  this  increase  must  be  due  to  an 
evolution  of  the  flame  temperature.  Taking  into  account  the  con¬ 
vective  velocity  of  the  flame  pockets,  which  is  roughly  5  m/s 
[18],  the  flame  pockets  have  a  lifetime  that  is  around  6  ms  longer 


Fig.  11.  Mean  and  RMS  of  temperature  in  flames  DJHC-I,  DJHC-V  and  DJHC-X  at 
z  =  90  mm.  The  fainter  upper  curves  indicate  T99.  the  99  percentile  of  temperature, 
to  indicate  measured  peak  temperatures. 


r  [mm] 

Fig.  12.  Mean  and  RMS  of  temperature  in  flames  DJHC-I, DJHC-V  and  DJHC-X  at 
z  =  120  mm.  The  fainter  upper  curves  indicate  T99,  the  99  percentile  of  temperature. 

at  z  =  120  mm  than  at  z  =  90  mm.  Apparently,  this  time  is  needed 
for  the  flame  pockets  to  reach  temperatures  close  to  equilibrium. 

The  DJHC-V  flame  does  not,  despite  its  higher  adiabatic  flame 
temperature,  show  higher  mean  or  peak  temperatures  at  z=  120 
than  DJHC-1.  Flame  structures  are  also  present  regularly  here 
(Pbl  rs  65%),  another  indication  that  chemistry  resides  far  from 
equilibrium. 

As  a  reference,  temperature  statistics  were  also  obtained  for  a 
conventional  lifted  flame  with  a  jet  Reynolds  number  of  5000,  at 
a  height  of  z  =  90  mm  (which  is  above  the  lift-off  height)  and  at 
the  radial  location  with  the  highest  mean  temperature 
(r=15mm).  The  pdf  of  temperature  of  this  flame  is  shown  to¬ 
gether  with  that  of  the  DJHC-I  flame  in  Fig.  13.  The  peak  tempera¬ 
ture  Tgg  of  the  lifted  flame  (indicated  with  the  right  dotted  line)  is 
2100  K.  When  rescaled  between  the  mixing  temperature  and  the 
adiabatic  flame  temperature,  this  is  comparable  to  the  Tgg  of 
DJHC-1  at  z=  120  mm.  The  temperature  distribution  of  the  lifted 
flame  is  much  wider,  containing  both  colder  samples  as  the  reac¬ 
tants  are  injected  at  room  temperature,  and  much  hotter  samples. 
Figure  14  compares  the  temperatures  at  z  =  90  mm  for  varying  jet 
Reynolds  numbers.  Only  minor  differences  can  be  spotted.  The 
mean  temperatures  vary  insignificantly,  the  Tgg  values  of  the 


T[K] 


Fig.  13.  Pdf  of  temperature  p(T)  (points  lie  on  the  center  of  the  bins)  for  the  DJHC-I 
flame  (Rejef  =  4100,  at  z  =  120mm,  r=14mm)  and  a  lifted  flame  (RcJf.;  =  5000  at 
z  =  90  mm,  r=  15  mm),  the  vertical  dashed  lines  denote  the  T99  values.  The  number 
of  measurements  is  951  (DJHC-I)  and  772  (lifted  flame). 
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Fig.  14.  Mean  and  RMS  of  temperature  for  three  jet  Reynolds  numbers  in  flame 
DJHC-I,  z  =  90  mm.  The  fainter  upper  curves  indicate  Tg9. 

lowest  jet  Reynolds  number  case  are  slightly  lower,  in  agreement 
with  its  higher  lift-off.  In  the  highest  jet  Reynolds  number  case, 
colder  temperatures  were  occasionally  measured  at  r=-17mm 
(2%  of  the  samples  were  below  1000  K),  originating  from  the  ambi¬ 
ent  air.  These  are  responsible  for  the  increase  in  the  RMS  values. 

The  differences  are  more  substantial  at  z=  120  mm.  Figure  15 
shows  that  the  mean  temperatures  as  well  as  the  peak  tempera¬ 
tures  are  substantially  lower  in  the  higher  Reynolds  number  case. 
The  mean  temperatures  in  this  case  are  even  lower  than  those  at 
z  =  90  mm,  which  can  only  be  explained  by  the  entrainment  of 
colder  coflow  air.  A  stronger  departure  from  chemical  equilibrium 
due  to  shorter  residence  time  of  flame  pockets,  and  possibly  the 
influence  of  larger  strains  could  form  an  additional  explanation 
for  the  fact  that  the  peak  temperatures  do  not  exceed  those  at 
z  =  90  mm  significantly,  as  they  do  in  the  lower  Reynolds  number 
cases.  As  the  influence  of  ambient  air  (witnessed  as  outliers  in 
the  pdfs  of  temperatures)  is  present  from  r=  12  mm  and  up  in 
the  high  jet  Reynolds  number  case  (in  the  other  cases  at 
r=20mm  and  up),  it  is  difficult  to  be  conclusive  regarding  this 
point. 

7.  Entrainment 

A  remarkable  feature  of  the  studied  flames  is  that  ignition 
events  are  first  observed  at  lower  heights  with  increasing  jet 


Fig.  15.  Mean  and  RMS  of  temperature  for  three  jet  Reynolds  numbers  in  flame 
DJHC-I,  z  =  120  mm.  The  fainter  upper  curves  indicate  again  T99. 


Reynolds  number.  An  explanation  for  this  phenomenon  will  be 
provided  in  this  section.  Of  interest  to  this  problem  are  the  proper¬ 
ties  in  the  mixing  layer,  located  radially  somewhere  between  the 
jet  centerline  and  the  pure  coflow  fluid  outside  the  jet.  Here,  mix¬ 
tures  of  fuel  and  coflow  fluid  with  the  optimum  mixture  fraction 
for  fast  ignition  are  formed  (in  general  not  equal  to  the  stoichiom¬ 
etric  mixture  fraction,  [26]).  The  total  mass  flow  in  the  jet  increases 
in  the  downstream  direction,  as  more  and  more  coflow  fluid  is 
drawn  into  the  jet.  Because  of  this  continuous  increase  in  mass 
flow,  the  coflow  fluid  that  is  entrained  by  the  jet  originates  from 
progressively  larger  radii  at  z  =  0  as  one  goes  downstream.  If  the 
properties  of  the  coflow  are  non-constant  along  the  radial  direc¬ 
tion,  the  properties  of  the  entrained  coflow  fluid  will  then  also 
change  in  the  downstream  direction.  For  the  DJHC  flames,  this 
means  that  the  fuel  stream  initially  mixes  with  the  colder  coflow 
fluid,  originating  from  small  radii  at  z  =  0.  As  more  coflow  fluid  is 
consumed  by  the  jet,  the  hotter  coflow  from  larger  radii  will  find 
its  way  into  the  mixing  layer  at  larger  heights.  The  temperature 
conditional  on  (the  most  reactive)  mixture  fraction  will  increase 
in  the  axial  direction  and  ignition  delay  times  will  become 
sufficiently  short  for  ignition  kernels  to  form  in  the  domain.  In  this 
analysis,  it  is  assumed  that  the  temperature  of  a  parcel  of  coflow 
fluid  as  it  is  entrained  is  related  to  the  temperature  of  that  fluid 
parcel  when  it  was  at  z  =  0.  This  is  justified  by  a  simple  scaling 
argument:  the  convective  heat  flux  is  much  stronger  than  the  dif¬ 
fusive  heat  flux  (uT  »  u"T"),  since  turbulent  fluctuations  are  rela¬ 
tively  small  in  the  coflow.  Because  the  length  scales  for  convection 
(length  of  streamline)  and  diffusion  (radius  of  the  burner  tube)  are 
of  similar  magnitude,  the  time  scale  for  turbulent  diffusion  is  much 
greater  than  the  timescale  for  convection.  As  a  consequence, 
temperatures  are  largely  conserved  along  streamlines  in  the 
coflow. 

Since  the  entrainment  flux  of  turbulent  jets  is  proportional  to 
the  centerline  velocity  (see  [27],  p.  35),  this  transport  of  hotter  co¬ 
flow  fluid  from  larger  radii  at  z  =  0  toward  the  jet  is  augmented 
when  the  jet  velocity  is  increased.  In  the  presence  of  a  positive  ra¬ 
dial  temperature  gradient  of  the  coflow  this  results  in  a  faster  rise 
of  temperatures  conditional  on  mixture  fraction  in  the  axial  direc¬ 
tion,  with  ignition  at  lower  axial  locations  as  a  result.  Note  that  in 
this  analysis,  two  important  assumptions  were  made.  First  of  all, 
details  of  the  mixing  process  itself  are  not  considered.  This  is  not 
to  say  that  these  complexities  (such  as  the  influence  of  scalar  gra¬ 
dients  on  the  mixing  and  ignition  process  [28])  are  not  present,  but 
rather  that  their  influence  is  overshadowed  by  the  effect  of  the 
temperature  variations.  Secondly,  ignition  delay  times  are  trans¬ 
formed  to  ignition  delay  lengths  using  a  constant  velocity,  inde¬ 
pendent  of  the  jet  velocity.  This  assumption  is  justified  by  the 
experimental  finding  that  the  convective  velocity  of  ignition  ker¬ 
nels  is  not  strongly  affected  by  the  jet  velocity,  but  is  mainly  deter¬ 
mined  by  the  coflow  velocity  [18]. 

The  transport  of  coflow  fluid  toward  the  mixing  region  in  the  jet 
can  be  studied  quantitatively  with  the  available  temperature  and 
velocity  data.  The  task  is  to  find  the  radial  position  r'  at  z  =  z0  from 
which  the  streamline  that  is  entrained  at  (r*,z*)  originates,  see 
Fig.  16.  By  definition  the  mass  flow  rate  through  any  boundary  that 
runs  from  a  point  p  to  any  point  q  on  a  streamline  S  is  constant, 
independent  of  the  choice  of  q.  If  the  value  of  the  streamfunction 
lF  is  set  to  zero  at  (r0,z0),  the  streamline  S('F)  in  Fig.  16  connects 
points  q'  and  q*  for  which  Mc  =  Me  =  ¥'. 

Now,  a  criterion  is  needed  to  define  the  interface  between  the 
turbulent  mixing  layer  and  the  coflow  stream.  This  criterion  is  pro¬ 
vided  by  the  turbulent  shear  stress  u'v’.  As  demonstrated  in  Fig.  10, 
the  value  of  this  shear  stress  at  z  =  60  mm  quickly  approaches  zero 
at  the  jet  edge,  enabling  the  use  of  a  stringent  threshold  level.  The 
edge  of  the  jet  is  defined  as  the  radius  r*  where  u’v’  is  less  than  1% 
of  its  peak  value. 
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Fig.  16.  Sketch  of  the  jet  boundary  (curved  thick  line)  and  the  streamline  S.  In 
absence  of  circulation  regions,  points  q'  and  q*  share  the  same  streamline  S,  when 
the  mass  flow  through  the  boundary  p-q*  (which  is  the  total  entrained  coflow  fluid 
by  the  jet  between  z0  and  z’)  equals  the  mass  flow  through  the  boundary  p  q’.  The 
thick  grey  line  denotes  the  boundary  between  the  jet  region,  where  jet  and  coflow 
fluid  mix,  and  the  coflow  region. 


The  value  of  the  entrainment  mass  flow,  Me,  can  be  obtained  by 
integration  along  the  jet  boundary,  as  was  done  by  Han  and  Mun- 
gal  [29],  This  requires  a  careful  determination  of  the  angle  of  the  jet 
boundary,  and  calculation  of  the  mass  flow  perpendicular  to  this 
boundary.  A  different  approach  is  taken  here,  using  a  control  vol¬ 
ume  that  runs  from  (r0,z0)  via  (r*,z0)  to  (r*,z*).  The  entrained  mass 
flow  between  z  =  z0  and  z  =  z*  is  the  sum  of  two  integrals,  one  in  ra¬ 
dial  and  one  in  axial  direction: 


Me(z*)  =  271 

/  pur  dr 

+  /  -pvrdz 

Jr0 

z„  K 

r*_ 

See  Fig.  16  for  a  definition  of  the  coordinates  r0,  r*.  z0  and  z*.  The 
mass  flows  in  Eq.  (5)  are  not  measured  directly.  Whether  and  how 
these  can  be  reconstructed  from  the  CARS  and  LDA  measurements, 
depends  on  the  nature  of  the  averages  of  these  measurements. 
Since  the  timing  of  CARS  measurements  is  controlled  externally, 
and  the  data  is  gathered  independently  of  processes  in  the  flame, 
the  temperature  measurements  can  be  considered  to  be  uncondi¬ 
tionally  sampled.  If  temperature  variations  within  the  measure¬ 
ment  volume  can  be  neglected,  the  arithmetic  mean  of  CARS 
measurements  therefore  corresponds  to  a  Reynolds  average.  For 
LDA  measurements  the  analysis  is  more  complicated.  In  Appendix 
A  it  is  shown  that  LDA  statistics  obtained  in  a  variable  density  flow 
with  application  of  transit  time  weighting  can  be  considered  to  be 
density  weighted.  Multiplication  of  this  density  weighted  mean 
velocity  by  the  mean  density  then  gives  the  mean  mass  flux,  as  re¬ 
quired  to  calculate  Eq.  (5),  since  by  definition  pit  =  pu  . 

This  point  is  not  of  critical  importance,  since  the  fluctuations  of 
density  are  small.  To  illustrate  this,  consider  the  maximum  value  of 
the  covariance  given  by  the  Cauchy- Schwarz  inequality: 

\p'v'\  <  (pl2vl2)1/2.  (6) 


The  density  variance  can  be  approximated  from  the  tempera 
ture  variance  using  a  linearised  relation  between  density  and  tern 
perature  around  the  mean  temperature: 

P  dT 


(7) 


If  the  measured  velocities  would  in  fact  be  Reynolds  averaged, 
the  uncertainty  (as  in  principle,  the  correlation  can  have  any  value 
between  one  and  minus  one)  is  typically  10%  or  smaller. 


Approximating  the  functions  in  the  integrals  of  Eq.  (5)  linearly 
between  the  measurement  points,  the  total  entrainment  can  be 
calculated.  The  totally  entrained  mass  flow  at  several  heights  is 
shown  in  Fig.  17.  The  entrainment  of  cases  DJHC-I  and  DJHC-X  with 
a  jet  mass  flow  of  16.1  nl/min  is  identical.  Furthermore,  the  higher 
velocity  jet  entrains  more  coflow  fluid,  as  can  be  expected  from  a 
standard  entrainment  scaling  [27].  The  data  shown  in  Fig.  17  can 
now  be  used  to  estimate  the  average  temperature  of  the  coflow  en¬ 
trained  by  the  jet  at  different  heights. 

The  values  of  T(r)  as  a  function  of  ¥' 0(r),  where  lF0  is  the  value 
of  Mc  at  z  =  3  mm,  are  plotted  in  Fig.  18.  The  coflow  temperature 
levels  off  at  around  f'o  =  0.6  x  10~3  kg/s,  corresponding  to 
rss  15  mm.  The  temperature  profile  is  much  flatter  against  F  than 
against  r,  because  the  majority  of  the  coflow  mass  flows  through 
the  center  and  outer  part  of  the  annulus. 

Comparing  now  Fig.  17  with  Fig.  18,  it  can  be  seen  that  at 
z  =  30  mm,  the  highest  jet  Reynolds  number  case  starts  entraining 
the  hottest  part  of  the  coflow  (at  around  1500  K),  whereas  this  hap¬ 
pens  at  around  20  mm  further  downstream  for  the  lowest  jet  Rey¬ 
nolds  number  case.  This  magnitude  is  in  agreement  with  that  of 
the  observed  decrease  in  the  height  where  ignition  kernels  are  first 
formed.  The  combination  of  the  dependence  of  the  entrainment 


Fig.  17.  The  total  entrainment  Me,  determined  by  numerical  evaluation  of  Eq.  (5),  as 
a  function  of  height  z,  for  four  different  cases.  The  errorbars  indicate  the  maximum 
contribution  of  the  density-velocity  covariance  calculated  with  Eq.  (6). 


Fig.  18.  Mean  temperature  as  a  function  of  the  streamfunction  ¥'  (set  to  zero  at 
r  =  3  mm),  at  z  =  3  mm. 
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fluxes  on  the  jet  velocity  with  the  positive  radial  temperature  gra¬ 
dient  of  the  coflow  therefore  forms  a  sound  explanation  for  the  ob¬ 
served  initial  decrease  in  lift-off  height. 

8.  Conclusions 

LDA,  CARS  and  OH-PLIF  measurements  were  performed  in  sev¬ 
eral  jet-in-hot-coflow  flames,  differing  in  jet  mass  flow  and  coflow 
temperature.  The  chemical  reactions  occurring  in  these  flames 
were  found  to  have  a  very  weak  impact  on  the  velocity  field,  as 
the  mean  velocities  and  turbulent  stresses  in  the  DJHC-I  flame 
are  comparable  to  those  of  the  DJHC-X  case,  that  has  a  lift-off 
height  exceeding  the  measurement  domain. 

In  the  temperature  field,  first  deviations  between  the  reacting 
and  non-reacting  cases  were  seen  at  z  =  90  mm.  These  differences 
are  not  very  pronounced  in  the  mean,  but  are  evident  both  in  the 
RMS  values  (nearly  doubling  when  reactions  start  to  occur,  case 
DJHC-I  vs.  DJHC-X),  and  in  the  99  percentile  of  temperature,  Tgg. 
The  peak  temperatures  in  the  DJHC-I,  Rejet  =  4100  case  at 
z  =  90  mm  are  far  away  from  the  adiabatic  peak  temperature 
whereas  those  at  120  mm  approach  this  temperature  much  closer. 
Considering  the  fact  that  flame  structures  are  already  present  most 
of  the  time  at  a  height  of  90  mm,  this  is  indicative  of  an  evolution 
in  temperature  of  flame  pockets  from  the  point  where  they  are 
newly  formed  to  where  they  have  evolved  longer  in  time  (a  differ¬ 
ence  in  time  of  around  6  ms).  In  other  flames  (DJHC-I  with  Re- 
jet  =  8800  and  the  colder  coflow  flame  DJHC-V  with  Rejet  =  4600), 
peak  temperatures  stay  remote  from  the  adiabatic  flame  tempera¬ 
tures,  also  at  z  =  120  mm. 

It  was  observed  with  both  high-speed  flame  luminescence  and 
OH-PLIF  measurements  that  the  location  of  first  chemical  reactions 
decreases  with  increasing  jet  velocity.  This  is  by  no  means  a  gen¬ 
eral  characteristic  of  jet-in-hot-coflow  flames,  but  is  related  to 
the  positive  radial  temperature  gradient  in  the  coflow  of  the 
DJHC-burner.  From  a  quantitative  analysis  of  the  entrainment  in 
the  near-nozzle  region,  it  was  concluded  that  the  hottest  part  of 
the  coflow  finds  its  way  into  the  jet  shear  layer  at  an  axial  location 
that  is  around  20  mm  lower  going  from  the  lowest  to  the  highest 
jet  values  of  the  jet  mass  flow.  This  distance  agrees  with  the  ob¬ 
served  decrease.  The  rate  with  which  hot  oxidizer  fluid  is  trans¬ 
ported  towards  the  edge  of  the  jet  where  mixing  takes  place  is 
thus  an  important  factor  in  the  flame  stabilisation  of  the  flames 
considered.  This  result  is  in  line  with  the  dominant  role  of  autoig¬ 
nition  in  the  stabilisation  of  jet-in-hot  coflow  flames  and  the 
known  sensitivity  of  autoignition  processes  to  temperature. 

In  the  Rejet  =  8800  case,  at  z  =  120  mm,  cold  outer  air  is  occa¬ 
sionally  entrained  in  the  location  where  the  flame  resides  as  seen 
by  outliers  in  the  temperature  pdfs,  highlighting  the  necessity  of 
using  a  larger  diameter  coflow  when  larger  jet  Reynolds  numbers 
or  higher  locations  are  to  be  studied. 


Appendix  A.  Density  weighted  nature  of  transit-time  corrected 
LDA  data 

Several  sources  of  statistical  biasing  in  LDA  data  have  been 
recognised,  most  notably  those  related  to  the  velocity  itself  [30]. 
In  flows  with  varying  density,  another  bias  will  arise.  Intuitively, 
one  can  see  that  a  fluid  parcel  with  low  fluid  density  will  have  a 
low  seeding  density  and  thus  a  lower  probability  of  detection  by 
the  LDA  system.  Therefore,  one  might  expect  the  velocity  statistics 
to  be  affected  by  these  density  fluctuations.  Goss  et  al.  [31  ]  showed 
that  in  their  combined  CARS/LDA  setup,  the  temperature  pdfs, 
conditional  on  the  presence  of  a  simultaneous  LDA  measurement 
(presence  of  a  seeding  particle),  were  very  similar  to  density 
weighted  CARS  measurements.  In  [32],  combined  LDA  and  CARS 


measurements  were  carried  out  to  investigate  velocity-tempera¬ 
ture  correlations,  using  the  LDA  measurements  to  compute  Favre 
averaged  quantities.  Nevertheless,  not  much  attention  has  been  gi¬ 
ven  in  literature  how  LDA  data  in  a  variable  density  environment 
should  be  interpreted. 

The  bias  in  LDA  measurements  occurs  because  the  data  rate  h 
(expected  number  of  seeding  particles  passing  through  the  mea¬ 
surement  volume  per  time)  is  a  function  of  the  value  of  the  instan¬ 
taneous  velocity  vector  (J,  and  possibly  of  other  properties  such  as 
the  density  p.  The  joint  pdf  of  the  measured  velocity  vector,  fu-LDA  is 
a  function  of  the  original  pdf  of  velocity  and  the  conditional  data 
rate: 

fu-LDA  =  fu  ,  (A.  1 ) 

with  the  averaging  operator  defined  as: 


Jv  fupdVdil/. 


(A.2) 


V  and  i //  form  the  sample  space  of  the  variables  U  and  p,  respec¬ 
tively.  The  instantaneous  data  rate  is  equal  to 

n  =  pcmV,  (A.3) 

with  V  the  volume  flux  and  cm  the  seeding  density,  [-/kg]  [33],  To 
eliminate  the  bias,  some  kind  of  correction  is  necessary.  This  is  real¬ 
ised  in  practice  by  applying  an  appropriate  weighting  factor  w  to 
the  data  points.  To  demonstrate  the  connection  between  applica¬ 
tion  of  weighting  factors  on  data  points  and  the  scaling  of  the  pdf, 
consider  the  variable  U  that  has  a  pdf  fu  and  is  weighted  with 
weight  factors  w  to  yield  the  weighted  pdf /U;w  : 


Here,  N  is  the  total  number  of  data  points,  and  M  the  number  of 
datapoints  between  14  —  LA u  and  14  +  j  Au.  Therefore,  if  the  condi¬ 
tional  expectation  of  the  weight  factor  is  inversely  proportional  to 
the  conditional  data  rate  (ii|V),  application  of  this  weight  factor 
on  Eq.  (A.l)  will  result  in  unbiased  statistics,  since  the  measured 
pdf  after  correction  is: 


f J±DAC  —  fu 


(h|V){w|V) 
(fiw)  ’ 


(A.5) 


Hoesel  and  Rodi  [33]  found  that,  in  a  constant  density  flow,  the 
transit  time  satisfies  this  criterion.  The  transit  time  is  a  random 
variable  for  a  given  value  of  V  (seeding  particles  pass  the  measure¬ 
ment  volume  at  arbitrary  locations),  with  expectation  : 


PM  = 


Vp 

V(V) 


Vppcm  p 
n(V)  x  ft(V)  ’ 


(A.6) 


with  Vp  the  probe  volume  (which  is  independent  of  the  flow  direc¬ 
tion,  see  the  Appendix  in  [33]  for  a  discussion).  The  LDA  data  re¬ 
ported  in  this  study  has  also  been  transit-time  weighted  to  avoid 
the  velocity  bias.  In  absence  of  density  fluctuations,  unbiased  esti¬ 
mates  for  mean  velocities  and  (co (variances  can  now  be  obtained 
by  weighting  each  measurement  with  its  individual  transit  time  tr,-. 

When  density  fluctuations  do  occur,  the  pdf  of  the  measured 
velocities  after  weighting  is  (using  (ii| V}  =  (p\V)cmV(V)  and  Eq. 
(A.5)): 


4  i 


(p\V)cmV(V)V-\V)  (p\V) 

,LDAL  =  JU - T32E3. - =  lu  - 


{p)  Cm 


(P) 


(A.7) 
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which  is  the  Favre  pdf.  Therefore,  if  transit  time  corrections  are  ap¬ 
plied  to  LDA  data  obtained  in  a  flow  with  uniform  seeding  density 
per  mass  (which  is  in  principle  not  influenced  by  thermal  expansion 
due  to  combustion),  density  weighted  quantities  will  be  obtained, 
i.e.: 


n  In 

I>fr'7  X><  = 

i=l  /  i=  1 

n  In 

^UiVitTi  / 
i=l  /  i= 1 


PU 

P 


puv 

~Y 


etc. 
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